Capillary filling in microchannels patterned by posts 
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We investigate the capillary filling of three dimensional micro-channels with sur- 
faces patterned by posts of square cross section. We show that pinning on the edges 
of the posts suppresses, and can halt, capillary filling. We stress the importance of 
the channel walls in controlling whether filling can occur. In particular for channels 
higher than the distance between adjacent posts, filling occurs for contact angles less 
than a threshold angle ~ 55°, independent of the height of the channel. 



I. INTRODUCTION 



Capillary filling, the ability of water to fill a hydrophilic channel, has been recognised since 
the pioneering work of Lucas and Washburn jl, 2], 3] nearly a century ago. However, inves- 
tigations of capillary filling in microchannels remain interesting due to modern applications 
of microfluid devices. Advances in lithographic techniques mean that it is becoming in- 
creasingly feasible to fabricate microchannels with well defined surface structures on micron 
ength scales. These have potential applications for chemical detection 4fl, as microreactors 
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5|, or to build entropic traps for DNA separation [6(. Our aim in this paper is to present a 
numerical investigation of how posts on the surface of a microchannel affect capillary filling. 
Our results are relevant to the use of electrowetting to control flow in microchannels and 
suggest ways to overcome the difficulties of filling structured microchannels. 

If a channel with hydrophilic walls comes into contact with a fluid reservoir it starts to fill 
as capillary forces pull the fluid into the channel. Balancing the capillary forces (27cos# a d) 
against the viscous drag of the entering fluid (12i]x(dx/dt)/H 2 ) gives an expression for the 
position of the advancing fluid in the channel function of time t 
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for a capillary of height H and infinite width. 7 is the gas-liquid surface tension, rj the 
fluid viscosity and 9^ the contact angle of the advancing front. This formula neglects 
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inertial and gravitational effects (good approximations once filling is established, and for 
channels of dimension smaller than the capillary length), assumes that the displaced fluid 
has zero viscosity, and neglects the slip length. These conditions are not always satisfied 
by simplified models used to investigation of capillary filling. However, excellent agreement 
between the theory and numerical results can be achieved by accounting for drag forces of 
the gas (12?7 eas (L — x) (dx / dt) / H 2 , where L is the length of the channel and 77 gas the viscosity 
of the gas) pj], or allowing for a slip length Q]. 
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FIG. 1: (Color online) (a) When a fluid interface reaches an edge it remain pinned for a range of 
angles 2tt — 29 eq — a. (b) A fluid front, advancing from the left, remains pinned (straight line) at 
the edges of two opposing ridges for all contact angles. 



When the surface of a microchannel is patterned with obstacles, such as posts or ridges, 
the capillary filling is, in general, suppressed because of pinning on the edges of the posts. 
The statics of the pinning can be understood by reference to the Gibbs' criterion 9|. This 
states that, when a fluid interface reaches a edge, it will remain pinned over the range of 
angles between the equilibrium contact angle on each of the surfaces bounding the edge as 
illustrated in Fig. [T^a). A striking consequence of the Gibbs' criterion is that if there are 



opposing ridges across the channel it does not fill [10j . Capillary forces pull the interface until 
it becomes a flat surface, pinned on the edges as shown in Fig. QJb). It cannot move down 
the sides of the post from this configuration, and hence will remain pinned. If the meniscus 
is moving, however, inertial effects may allow it to overcome the pinning as it reaches the 
ridge. Then, as the interface moves down the channel it slows, and will eventually pin at a 
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subsequent ridge. We also note that for nanoscale roughness thermal fluctuations allow the 
interface to depin and advance 

If the ridges across the channel are replaced by separated posts capillary filling becomes 
possible for sufficiently low values of the contact angle. The aim of this paper is to investigate 
how the depinning behaviour depends on the channel geometry and the contact angle, and 
to explore the depinning mechanism in some detail. We emphasise the importance of the 
channel walls in controlling the depinning of the interface and we identify the way in which 
the interface depins for different channel geometries. Estimates are given for the contact 
angles at which depinning occurs. 

In Sec. [Til we define the model we use, summarising the equilibrium properties and the 
equations of motion. We then describe the geometry of the channel and list the simulation 
parameters of the model. In Sec. II III we present our results. For two typical channels, 
using values of 8 eq which do not pin, we compare the filling rate to the similar case of a 
smooth channel (without obstacles), showing how posts on the channel surface slow the 
filling. We then concentrate on the filling/pinned transition. The case in which the posts 
span the channel and that in which the posts do not meet across the channel are investigated 
in Sees. IIII Al and IIII Bl respectively. We find that, for high channels and long posts, the 
pinning-depinning threshold angle 6^ ~ 55°. The depinning is driven by the walls and is 
independent of the channel geometry. For narrow channels or short posts, filling is possible 
for higher contact angles. In Sec. IIII CI we discuss the effect of inertia on the determination 
of the threshold contact angle 9^. Finally Sec. [TV] presents our conclusions. 



II. THE MODEL 



As we are considering micron length scales it is appropriate to describe the system using 
a mesoscale modeling approach. We choose to use a diffuse interface model, solved using 
a lattice Boltzmann algorithm Jj], Q] , which has proven to be a useful tool to model 
the dynamics of fluids with moving interfaces. Of particular relevance here Ref. [15| demon- 
strates how this approach can be used to describe capillary filling in smooth microchannels. 
We now give details of the model, the channel geometry and the simulations parameters used 
in this paper. Details of the implementation of the lattice Boltzmann method of solving the 
equations can be found in [l^], and are not repeated here. 
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A. Equations of motion 



We consider a binary fluid with components A and B, say, described by the free energy 
functional 



—n\ogn + -(d a <p) --(f) +-c 



+ 



h- 



(2) 

where n is the local total density of the A and B components (n = + (f) is the 
order parameter <p = ha — ns and c is the lattice velocity c = 5x/St, where Sx is the lattice 
spacing and 5t is the simulation time step. The first integral in Eq. ([2]), taken over the total 
volume Q, controls the bulk properties of the system. The terms in cf) give coexistence of 
phases with = ±1. The energy cost of an interface between the two phases is modeled 
by the derivative term, with k related to the surface tension. The term in n controls the 
compressibility of the fluid. 

The integral over the solid-liquid interface dQ in Eq. ([2]) accounts for the wetting prop- 
erties of the solid surfaces, h is related to the equilibrium contact angle 8 eq by fl6 | 
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with sign(a;) = 1 if x > and sign(x) = — 1 otherwise. 

The hydrodynamics of the fluid is described by the the Navier-Stokes equations for the 
density p and the velocity field v together with a convection-diffusive equation for the binary 
order parameter <p 

d t p + V- (pv) = , (4) 

d t {pvp)+d a {pv a vp) = -d a [P aP + r}{dpv a + d a vp)) , (5) 

d t (j) + V • (0v) = MV 2 /i . (6) 

In Eq. (J5j) rj is the viscosity of the fluid and in Eq. |6} M is a mobility coefficient. The 
pressure tensor P a p and the chemical potential \i which appear in Eqs. (jSJ) and (jSJ), which 
describe the equilibrium properties of the fluid, follow from the free energy (J5J) as 



a/3 



On 



5, 



a/3 



8^ 

V = TT 



+ 5, 



a/3 



5^ 5^ 

3a 



3 n + T 4 - -0 2 - -{d T 4>Y ~ K<f>d T d T 



a0 3 2 — K,d T d T 



(7) 
(8) 



5 



We have chosen to use a two- component binary fluid as a model system to simulate 
capillary filling. This is because modeling capillary filling correctly using a liquid-gas diffuse 
interface model is computationally demanding because of unphysical motion of the interface 
due to evaporation-condensation effects |7j. However our results are equally applicable to a 
physical system where a liquid displaces a gas as the important physical parameters are the 
viscosities, not the densities, of the fluid components. Therefore we shall use the natural 
terminology 'liquid' and 'gas' for the displacing and displaced fluid from now on. 

B. Simulation geometry 




X 



(a) (b) 

FIG. 2: (Color online) Structure and parameters of the channel. The fluid flows from left to right 
along the x-axis. H is the height of the channel, H e ^ the distance between the tops of the posts, 
D is the dimension of the square posts and L is the distance between the posts. Fig. [2ja) is a cross 
section in the x—z plane, and Fig. EJ^b) is a three-dimensional view of the channel geometry. 

Fig. [2] reports the channel geometry which we use in this paper. Two reservoirs, of liquid 
and gas, which are in contact to equalise the pressure, are connected to a channel, running 
along x, with walls decorated by equispaced rectangular posts. The relevant geometric 
parameters are the channel height H (measured from wall to wall), the distance between 
the top of two posts on opposing walls H e s, the distance between two obstacles across the 
channel L, and the cross section of the posts which we choose to be square with side of 
length D. 
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C. Simulation parameters 

All the quantities reported in this paper are expressed in units of Sx, the lattice spacing, 
and St, the time step and hence c = 1. We simulate channels with one, two or three rows 
of posts corresponding to lengths (in the x direction, see Fig. [2]) L x from 80 to 200 lattice 
spacing. In the y direction we employ periodic boundary condition and sizes L y from 60 to 
170, while H = L z spans from 30 to 100. For the reservoirs we use the same L y res = L y 
while in the z direction we use periodic boundary conditions with L z>Tes two or three times 
L z . The x- length of the reservoir L XtTes is approximately half the channel length L x . For the 
free energy (j2J), the bulk phases (0 = 1 and <fi = —1) are interpolated by an interface with 



a profile which is well approximated by = tanh(:r/\/2£) (with £ = ^Jk/a), and with a 
surface tension equal to 7 = ^JSka/9. Here we use a = 0.04 and k = 0.02. These values give 
an interface width of order four lattice Boltzmann nodes which is much smaller than the 
typical size of the posts used. The viscosity rj and the mobility coefficient M appearing in the 
hydrodynamic equations (jSl E]) are related to the relaxation time in the lattice Boltzmann 



algorithm 



12j. We use M = 0.5. For the gas viscosity ?7 gas = 0.033, while for the liquid 



(unless specified) r]\ iq = 0.83. Both the gas and liquid densities are set to one. 



III. RESULTS 



In a typical simulation we start with a configuration in which the liquid-gas interface 
has advanced to the end of the first row of posts as indicated by the broken line in Fig. 
[2j If the contact angle is not too large, the interface moves through the channel, driven by 
the capillary force. If the filling fluid overcomes the first row of obstacles, we have verified 
it depins from the second row too. We have also used starting configurations with empty 
channels, finding no differences in the filling/pinned phase diagram. 

Fig. [3] shows the position of the advancing front (in the middle and near the walls of 
the channel) as a function of time for two typical geometries and for contact angles which 
do not pin. It is immediately apparent that the flow profile is very different to that of a 
smooth channel H e $ = H . Three regimes are present. When the advancing fluid reaches 
the beginning of the obstacles (dotted lines in Fig. [3]) it accelerates because of the increase 
in the capillary force due to more wettable surface provided by the obstacles. When the 
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FIG. 3: (Color online) Position x of the advancing front for (a) D = 20, L = 40, H eS = 0, H = 40, 
d eq = 55°; (b) D = 20, L = 40, H cS = 30, H = 50, 6 cq = 60° as a function of time t. The position 
of the advancing front near the wall (full line) and in the middle of the channel (broken line) are 
compared. We compare the filling rate for the same channel without posts (-ff e ff = H). 

front reaches the end of the obstacles it remains pinned for a certain time during which it is 
almost at rest. When finally it depins it restarts filling the channel with a Lucas Washb urn- 
like law (pQ), but with a reduced velocity (the "decelerated" regime compared to the smooth 
channel in Fig. [3]). This happens because the drag force is now larger due to the presence 
of the obstacles within the displacing viscous fluid. We observe that the "accelerated" and 
"decelerated" regimes are more obvious for if eff = 0, because the post surface is bigger than 
for H cS = 30. 

For the parameters considered in Fig. [3] the front is finally able to depin from each of the 
obstacles to move down the channel, but for lower contact angles it remains pinned. We 
next look in more detail at the pathways for depinning. 

A. H eS = 

We first consider geometries in which the posts span the channel (H e g = in Fig. [2]). Fig. 
H] summarises results which distinguish the cases where the interface is pinned on the posts 
from those where it can advance along the channel, for a range of geometric parameters (H, 
L and D, denned in Fig. [2]) and equilibrium contact angles 9 eq . 
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FIG. 4: (Color online) Filling/pinning transition for the H e ^ = geometry (posts that span the 
channel, see Fig. [2]). Open symbols are simulations in which the front is pinned while full symbols 
correspond to simulations in which the channel is filled. For H/L > 1 the filling-pinning transition 
happens for a threshold angle ~ 55° while for H/L < 1 the transition is possible for a range of 
equilibrium contact angles which increases with decreasing H/L. 

At high 6 eq the front remains pinned. However, at lower 9 eq the meniscus can overcome 
Gibb's pinning and the channel fills. This is due to the presence of walls bounding the 
channel. By advancing along the walls the meniscus is able to reach the angle it needs to 
move across the face of the posts. This is illustrated in Fig. [5] for different aspect ratios of 
the channel. 

Two different regimes are apparent in Fig. [2J For H/ L > 1 the boundary between contact 
angles that allow filling and those that do not is independent of H/L occurring at an angle, 
that we shall denote 9^, ~ 55°. For H/L < 1, however, the transition angle is not constant 
but increases with decreasing H/L. 

Fig. [5] compares the way in which the front depins from the posts in each of the two 
regimes. For H/L > 1 (Fig. [5] column (a)) the walls act independently. The menisci from 
two neighboring gaps first meet at the walls and then the advancing front covers the posts, 
moving from the walls towards the centre of the channel. Hence the contact angle below 
which depinning proceeds, 8^, is independent of H/L. 
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FIG. 5: (Color online) Depinning path ways for H eS = 0: (a) H/L > 1 (D = 20, L = 40, H = 80, 
6i cq = 50°) and (b) H/L < 1 (D = 20, L = 150, H = 50, 8 cq = 70°). The first three rows are three 
dimensional views of the advancing front (aspect ratios not to scale) and the fourth row shows the 
position of the advancing front on the face of the posts as a function of time during depinning. 
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FIG. 6: (Color online) Position of the depinning front near the walls (at z = 2), as a function of 
time (for £ = 1,2 and 3 • 10 4 lattice Boltzmann time steps) for 9 eq = 50° and several H so that 
H/L > 1 with L = 40. 

Further evidence is provided for this by Fig. [6] which shows how the meniscus advances 
along the wall for different values of H and at fixed L. The profiles are nearly independent 
of H as long as the front remains pinned. Without obstacles, Lucas- Washburn's law predicts 
a velocity proportional to H (Eq. ([I])), so that the profiles in Fig. would be well separated 
in a smooth channel. 

We stress the importance of the walls even in the limit H — > oo. Even for a very high 
channel depinning will occur at 9^, and will proceed as in Fig. [5] column (a) (except that it 
will take more time to wet the posts). However, without any walls the advancing fronts will 
remain pinned at the obstacles and flat as in Fig. [U(b). 

For H/L < 1, however, depinning occurs for contact angles greater than 6^. This is 
because the walls are sufficiently close that the interface moves in a concerted way across 
the channel. Therefore, once the menisci have advanced along the surfaces sufficiently far 
for depinning to occur, the interface depins along all of a post at the same time, and the 
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posts are wet from the sides. As H/L decreases the two surfaces more easily deform the 
interface and hence depinning can take place at a higher contact angle. This agrees with 
analytic results showing that, in the H/L — > limit, the interface depins for all hydrophilic 



contact angles 



171 ]. and can be understood using a free energy argument. The free energy 



gain in deforming the advancing interface scales like the area of the gas-liquid interfaces, 
and goes to zero as H for H — > 0. However the loss in free energy due to the wetting of 
the walls remains constant for H — ► 0. As a consequence the interface will advance for any 
9 cq < 90° for sufficiently small H. 

We expect the phase diagram to depend only on H/L and to be independent of D as, 
once the front has started to depin, the interface will continue moving across the post until 
it covers it. Indeed, using L = 40 we repeated simulations for D = 30 (H = 60, 80) and 
D = 40 (H = 60, 80) finding, again, a depinning transition compatible with On the 
other hand we observed a weak dependence of on L. We repeated simulations for L = 30, 
D = 30 and H/L = 20, 30, 50, 80, 100 and 120. 0{£ was again independent of H/L for 
H/L > 1, but placed between 50° and 55°, while for L = 40 (Fig. E> 9 eq = 55° fills the 
channel. In Sec. IIII CI we show the dependence of 6^ on L is related to inertial effects. 



B. H cS > 

In this section we describe the behaviour for the more general case H c g > when the posts 
do not reach all the way across the channel. We will consider the H/ L > 1 geometry where 
the interface depins from posts with H c g = by a transition at 9^, driven by depinning 
initiated at the walls. For H e ^ > filling the channel should be easier. 

Fig. [7] shows pinned configurations of the meniscus near the channel wall for different 
values of H e s, keeping the other geometric parameters fixed, and for 9 cq = 60°. For H c s = 
this corresponds to a pinned configuration. The effect of increasing H e g is to increase the 
contact angle at the edge of the post. When this angle exceeds the equilibrium contact angle 
(# eq = 60° in Fig. [7j) then, as predicted by the Gibbs' criterion, the front depins. For the 
parameters of Fig. [7] this happens for H e g = 40 (not reported in the figure). 

For H e g < H/2, the angle the meniscus makes with the edge of the post does not change 
significantly. For this reason, at low values of H e s the threshold equilibrium contact angle 
is almost the same as in the H e g = case (i.e. ~ 9^). However the meniscus advances 



12 





H e«= 2 
- H e«= 30 
■ H e«= 34 







10 20 30 40 50 

y [simulation units] 



60 



FIG. 7: (Color online) Position of the pinned meniscus near the walls (z = 2) for H = 60, D = 20, 
L = 40, 6> eq = 60° and four values of H c r. Increasing H c r the contact angle at the edge of the post 
is increased until the depinning threshold is reached (for H e R = 40). 



significantly from its position at H c q = when H e g becomes comparable to H. This can be 
explained by considering the shape of the pinned interface. If if e ff = and H/L > 1, the 
pinned front is almost flat in the middle of the channel with significant deviations only near 
the walls. If H e r is small enough that the gap does not overlap with the wall regions, then 
it will not perturb the interface compared to the H c s = configuration. On the other hand 
for large H e g, the end of the posts will lie within the wall region and there will be significant 
deformation of the pinned front, compared to the H e R = case, and a threshold equilibrium 
contact angle different to 9^. 

For several channels (D = 20, L = 40, H = 40, 50, 60, 70, 80) we have verified that 
at low H e s the front remains pinned at 9 eq = 60° > 9^, as in the H eS = case, while at 
high enough H e ff, 9 eq = 60° can fill the channel. We were never able to observe filling for 
= 65°, but we do not exclude that at high enough H eS (which would require a larger 
simulation box to properly resolve the small height of the posts) this can happen. 
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FIG. 8: (Color online) Full symbols denote values of H c g (with D = 20, L = 40 and H = 
40,50,60,70,80) for which the channel will fill at an equilibrium contact angle 9 eq = 60°. Empty 
circles, correspond to interfaces that remain pinned at 9 eq = 60° as in the case H e R = 0. 

Results summarising the threshold in H e g above which 8 eq = 60° fills are reported in 
Fig. [HI Having fixed D and L, we use H e ^/H and H/L as control parameters. At low value 
of H/L the filling for 8 cq = 60° occurs for values of H e g/H greater than 0.5 (see full circles). 
Increasing H/L higher values of H c r/H are needed to achieve filling at 60° because the 
length of flat interface in the middle of the channel increases with increasing H . 

At small values of H e g (for which 6^ = 60° does not fill) the depinning route is the same 
as for H e ff = 0. This is shown in the first column of Fig. [91 As the meniscus advances near 
the walls of the channel the interface between the post remains almost completely flat. The 
posts are wet, as before, from the channel walls towards the centre. This is not the case for 
larger value of H e R (full symbols in Fig. [8]). In Fig. [9] column (b) the advancing front does 
not remain pinned in the middle of the channel but the post are wet from their ends and 
sides towards the walls. 
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C. Inertial effects 

In performing the simulations we found that inertial effects made it difficult to determine 
the exact contact angle at which depinning occurs. If the front reaches the free energy 
minimum which corresponds to pinning with a residual kinetic energy it can overshoot, and 
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FIG. 9: (Color online) Depinning mechanism for (a) H = 80, H e g = 30, L = 40 and # e q = 55°; 
(b) = 60, H e f[ = 40, L = 40 and # e q = 60°. In the first row we report three dimensional 
interface profiles during depinning, while in the second row the position of the advancing front 
along a section taken through the centre of a post at different times. In (b) the face of the post is 
wet both from its ends and sides. 



hence depin, and it is not possible to entirely eliminate this effect without prohibitively long 
simulations. Indeed in a physical system the front will approach the posts with a finite 
velocity and whether it will pin will be a balance between 9 eq and the extent to which the 
front has been slowed by the viscous drag in the channel. 

To demonstrate the effects of inertia close to the depinning transition we investigated 
filling a channel with D = 20, L = 40, H = 60,if e fr = 0, # e q = 55°, from the starting 
configuration shown in Fig.[2J for three different liquid viscosities. Fig. [TD] shows the position 
of the advancing front as a function of time for each case. At early times lower viscosities give 



15 
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FIG. 10: Position of the maximum of the advancing front near the wall {z = 2) as a function of 
time for D = 20, L = 40, 9 eq = 55° and three different values of the viscosity rj. 

slightly higher speeds. However, the higher viscosity fluids gain more energy from the walls 
and have sufficient inertia to move past the pinning position for this contact angle whereas, 
at least on the timescale of the simulation, the lowest viscosity fluid remains pinned. As 
anticipated in Sec. IIII A\ inertial effects are also responsible for the dependence of on L, 
which controls the amount of water advancing between two posts (or equivalently the scale 
of the system) . 

We tried to better determine the position of 9^ repeating simulations in which 9 eq was 
gradually decreased (by 0.5° each 10 5 time steps) to better reproduce a quasi-static relaxation 
of the interface. For a channel with L = 30 we obtained 9^ = 52.5° ± 0.5°. However this 
estimate is very difficult because of the flatness of the free energy profile and the very slow 
interface velocities. 

We stress that, for angles near depinning, inertial effects depend primarily on the way 
in which the interface is pulled beyond the end of the posts by the surface, but not on its 
initial position within the microchannel. Hence we expect that this is not just an artifact of 
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the simulations, but that similar inertial effects will occur in experimental systems. 

IV. DISCUSSION 

As microfabrication techniques become standard it is becoming possible to design mi- 
crochannels with complicated internal geometries that may prove useful in controlling fluid 
behaviour. As a step towards understanding how fluids move in such channels we have in- 
vestigated capillary filling in microchannels patterned by regularly spaced square posts. A 
consequence of the Gibbs' criterion is that ridges that face each other across a channel will 
always pin a slowly moving interface. We show that, if the ridges are replaced by posts, the 
interface is able to depin for sufficiently small contact angles. This is because the meniscus 
can advance along the surfaces of the channel between the posts, thus allowing depinning 
to occur. 

For posts which span the channel and for a ratio of channel height to distance between 
the posts H/L > 1 the depinning threshold 9^ is independent of H/L because the two 
surfaces of the channel act independently. 9^ lies between 55° and 60° and the posts are 
wet from the surfaces towards the centre of the channel. For Hj L < 1 the threshold contact 
angle increases with decreasing H/L as the surfaces act cooperatively to reduce interface 
curvature across the channel. Here the posts primarily wet from their sides to their centres. 

In the general case in which the posts on opposing sides of the channel are separated 
by a distance H e g two regimes are present for H/L > 1. At low H e g the filling/pinned 
transition is similar to that for H e g = 0, with a threshold value of the equilibrium contact 
angle around 9^ and posts which wet from the wall to the centre of the channel. For high 
enough i? e ff, however, the threshold equilibrium contact angle for depinning increases, and 
the posts can also wet from the centre of the channel towards the walls during depinning. 



These results are in agreement with [lOj , where channels with fixed H and if e ff and several L 
were considered. In particular in the range H/L > 1, a contact angle threshold compatible 
with 9^ was observed. On increasing the value of L (i.e. exploring the H/L < 1 regime), 
filling at higher 9 eq was found, in agreement with our H e r = results. 

We have concentrated mainly on the quasi-static situation where inertia is neglected and 
therefore our threshold values are relevant to a very slowly moving interface. We have, 
however, shown that close to the threshold even tiny interface velocities can aid depinning. 
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It will now be interesting to investigate the more general case of a moving interface and 
assess the extent to which dissipation at the posts can slow and eventually pin the interface 
for a range of contact angles. 

The contact angle of a fluid within a microchannel can rather easily be varied by applying 
an electrowetting potential. This opens the possibility of controlling the fluid motion by 



switching 9 cq in and out of the pinning regime [18f|. This is of particular interest in the 
large H/L regime because the threshold equilibrium contact angle is well approximated by 
9^, independent of the channel geometry. 
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